Freezing and melting of 3D complex plasma structures under microgravity conditions 

driven by neutral gas pressure manipulation 
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Freezing and melting of large three-dimensional complex plasmas under microgravity conditions is 
investigated. The neutral gas pressure is used as a control parameter to trigger the phase changes: 
Complex plasma freezes (melts) by decreasing (increasing) the pressure. Evolution of complex 
plasma structural properties upon pressure variation is studied. Theoretical estimates allow us to 
identify main factors responsible for the observed behavior. 
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Complex (dusty) plasmas - systems consisting of 
highly charged micron-size particles in a neutralizing 
plasma background - exhibit an extremely rich variety 
of interesting phenomena [1]. Amongst these the transi- 
tions between fluid and solid phases (freezing and melt- 
ing) are of particular interest [1-4]. This is largely a 
consequence of the fact that high temporal and spatial 
resolution allows us to investigate these phase changes 
along with various related phenomena at the individual 
particle level [H-Q- 

In this Letter we report on experimental studies of 
the sohd-liquid phase changes in large 3D complex plas- 
mas under microgravity conditions driven by manipulat- 
ing neutral gas pressure. It is observed that the system 
of charged particles exhibits freezing (melting) upon de- 
creasing (increasing) the pressure, in contrast to the sit- 
uation in ground-based experiments where plasma crys- 
tals melt upon reducing the pressure. This can illustrate 
important differences between generic (e.g. similar to 
colloidal suspensions) and plasma-specific mechanisms of 
phase transitions in complex plasmas. 

Experiments are performed in the PK-3 Plus labo- 
ratory operating onboard the International Space Sta- 
tion [8|. The heart of this laboratory is a parallel-plate 
radio- frequency (rf) discharge operating at a frequency 
of 13.56 MHz. Discharge can operate in argon, neon or 
their mixture in a wide range of pressures, rf-amplitudes 
and rf-powers. Complex plasmas are formed by inject- 
ing monodisperse micron-size particles into the discharge. 
The particles are illuminated by a diode laser and the 
scattered light is recorded by CCD cameras. For a com- 
prehensive review of the PK-3 Plus project see Ref. [8]. 

The experiments described here are carried out in ar- 
gon at a low rf-power (~ 0.5 W). We use two differ- 
ent sorts of particles in the two distinct experimental 
runs: Si02 spheres with a diameter 2a = 1.55 jam and 
Melamine- Formaldehyde spheres with 2a = 2.55 /im. 
The experimental procedure, identical in these two runs, 
is as follows: When the particles form a stable cloud in 



the bulk plasma, the solenoid valve to the vacuum pump 
is opened, which results in a slow decrease of the gas pres- 
sure p. Then, the valve is closed and p slowly increases 
due to the gas streaming in. (Neutral flow has negligi- 
ble direct effect on the particles). During the pressure 
manipulation {c^ 6 minutes in total), the structure of the 
particle cloud is observed. The observations cover the 
pressure range from 15 Pa, down to the lowest pres- 
sure of 11 Pa and then up to 21 Pa (see Fig. [TJi). 
In order to get three-dimensional particle coordinates, 30 
scans are performed. Scanning is implemented by simul- 
taneously moving laser and cameras in the direction per- 
pendicular to the field of view with the velocity 0.6 mm/s. 
Each scan takes 8 s, resulting in the scanning depth 
of 4.8 mm; the interval between consecutive scans is 
4 s. The particle positions are then identified by tomo- 
graphic reconstruction of the 3D-pictures taken with the 
high resolution camera, which observes a region of about 
8x6 mm^ slightly above the discharge center. 

Let us first analyze the global reaction of the parti- 
cle cloud on pressure manipulation. Figure [1^ shows 
the positions of the upper and lower cloud boundaries 
and the cloud thickness as a function of the scan num- 
ber (time). It is evident that the position of the upper 
boundary is strongly correlated with pressure: It moves 
downwards (upwards) with the decrease (increase) of p. 
This has a clear physical explanation. Particles cannot 
penetrate in the region of strong electric field (sheath) 
established near the upper electrode. The position of the 
upper cloud boundary is thus set by the sheath edge. The 
sheath thickness is roughly proportional to the electron 
Debye radius Adc which exhibits the following approxi- 
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upon a decrease (increase) of p, the particles are pushed 
farther (closer) to the electrode, in full agreement with 
the observations. The lower cloud boundary, associated 
with the presence of the particle-free region (void) in the 
central area of the discharge [sl, llO|] , shows less system- 
atic behavior. Its position does not change when p de- 
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FIG. 1: (Color) (a) Positions of the upper and lower bound- 
aries of the particle cloud and its thickness in the vertical di- 
rection vs. the scan number for the experimental run with 
small particles. The corresponding values of pressure are 
shown for both runs [red dashed (blue solid) curves for small 
(big) particles]. Left inset shows side view of the particle 
cloud (inverted colors), right inset corresponds to the FoV 
of the high resolution camera [particles are color-coded to 
see solid-like (red) and liquid- like (blue) domains]. Rectangle 
marks the part of the cloud used for the detailed structural 
analysis (rectangular box 7 x 0.7 x 4.5 mm^); (b) Interparticle 
distance vs. scan number in two experimental runs. Insets 
show the dependence A(p) demonstrating some hysteresis. 

creases, but then moves upwards monotonously when p 
increases. However, the displacement amplitude is rela- 
tively small. As a result, the thickness of the cloud ex- 
hibits pronounced decrease (increase) when p decreases 
(increases), see Fig. [TJi. Thus, the particle component 
becomes compressed by reducing the pressure. The re- 
sulting dependence of the mean interparticle distance (in 
the part of the particle cloud subject to detailed anal- 
ysis) on the scan number /pressure is shown in Fig. [TJd. 
The interparticle distance A is clearly correlated with 
p, although some hysteresis (more pronounced for small 
particles) is evident from insets in Fig. [TJd. 

For the detailed analysis of the structural properties 
during the pressure manipulation, a part of the cloud 
sketched in Fig. [T] has been chosen. It contains of the or- 
der of 10^ particles. To determine the local order of par- 
ticles we use the bond order parameter method 11 ll . In 
this method, the local rotational invariants g^, Wi |lll4l3l| 
for each particle are calculated and compared with those 
for ideal lattice types. Here, we are specifically interested 
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FIG. 2: (Color) Variation of the structural properties with 
pressure as reflected by particle distributions on the plane 
of rotational invariants (^4,^5) for the system of small (left 
panel) and big (right panel) particles. The rotational invari- 
ants for perfect hep and fee lattices and liquid-like domain 
[sketched in (a)] are also indicated. For discussion see text. 

in identifying face-centered cubic (fee), body-centered cu- 
bic (bcc), and hexagonal close-packed (hep) lattice types 
and, therefore use the invariants ^4, qq^ wq calculated 
using 12 and 8 nearest neighbors. Figure [2] shows repre- 
sentative examples of particle distributions on the plane 
(^4, Qe)' Initially, both systems of smal (Fig. [2^) and 
big (Fig. [2li) particles reveal weakly ordered liquid-like 
phase. Upon a decrease in p, the particles tend to form 
more ordered structures, and eventually both systems 
freeze. The distributions corresponding to the minimum 
pressure are shown in Figs.[2l3, e. Clear crystalline struc- 
tures which are dominated by the hep and fee lattices 
are observed. Subsequent increase in pressure destroys 
the particle ordering. Figures [2] c,f demonstrate the final 
states of the systems. 

Scan by scan analysis of the complex plasma struc- 
tural composition reveals the following main properties: 
(i) A decrease (increase) in p enhances (suppresses) or- 
dering of the particles; (ii) Maximum number of particles 
in the crystalline state corresponds to the minimal pres- 
sure; (iii) Crystalline phase is mostly composed of hcp- 
and fcc-like particles with only a small portion of bcc- 



FIG. 3: (Color online) Distributions of small particles over 
the rotational invariant wq for different scan numbers (indi- 
cated in the figure). The corresponding cumulants Weix) are 
plotted in the upper part of the figure. The position of the 
half-height of We{x) is used as a melting indicator (see text). 



like clusters; (iv) The system of small particles exhibits 
melting with increasing whilst the system of big parti- 
cles approaches the melting transition, but does not fully 
melt; (v) Premelting stage (disappearance of fee- and 
bcc-like particles [12|) is observed for small particles, but 
is less pronounced in the system of big particles. 

It has been shown recently that cumulative distri- 
bution function Wq{x) = J_^n{wQ)dwQ is very sen- 
sitive to phase state of various systems 0, Here 
n{wG) is the distribution, normalized to unity, of par- 
ticles over the rotational invariant wq. Figure [3] shows 
n{wQ) and the corresponding W6{wg) for four different 
scans corresponding to different states of the system of 
small particles. An appropriate melting indicator can 



be defined as = Wq^/wq^"^^ where Wq^ is the posi- 
tion of the half-height of Weiwe) [Weiw^^) = 1/2] and 



,,fcc 



-1.3161 X 10- 



is the value of Wq^ for the fee lat- 



tice. Melting occurs when ^1-3. There also exists a 
freezing criterion [l5| based on the properties of the pair 
correlation function g{r). It states that near freezing, the 

ratio R = ^min/^max^ ^j^g^g ^max ^^min) -g ^^i^^ 

g{r) at the first maximum (first nonzero minimum), is a 
constant, R o:^ 0.2. This rule describes fairly well freez- 
ing of the Lennard- Jones fluid, but is not universal: R 
can vary considerably for repulsive potentials of differ- 
ent softness [16]. However, softness of the interactions 
in the considered complex plasmas is expected to fall in 
the range, where R is nearly constant. Fig. shows 
melting and freezing indicators P^ and R for different 
scans. Both criteria clearly indicate that the structures 
are mostly ordered near the lowest pressure and identify 
freezing and melting at similar conditions [except system 
of big particles, which does not melt (is near melting) 
according to the value of R (P^)]. 

We interpret the observed fluid-solid phase changes to 
the variation in the electrical repulsion between highly 
charged particles. Manipulating the gas pressure exper- 
imentally, changes various complex plasma parameters 
and modifies the strength of the repulsion. When the 
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FIG. 4: (Color online) (a) Values of the melting (freezing) 
indicator P^ {R~^) vs. the scan number. Solid curves cor- 
respond to P^, dashed curves to . System of small (big) 
particles is denoted by blue (red) color, (b) Estimated inter- 
particle interaction energy for each scan. 



electrical coupling reach (or drops below) certain level, 
freezing (or melting) occurs. To verify this scenario we 
need to estimate relevant plasma parameters. 

We use the results from siglo-2D simulations [8] to 
estimate plasma parameters in the absence of parti- 
cles. In the considered regime {p ^ 10 — 25 Pa, rf- 
amplitude ~ 15 V) the central plasma density is linear 
on p and can be, with a reasonable accuracy, described 
as no — (1.20 + O.llp) x 10^, where no is in cm~^ and 
p in Pa. The electron temperature exhibits almost no 
dependence on p, Tg 3.8 eV. Ions and neutrals are 
at room temperature Ti^n ^ 0.03 eV. When particles 
are injected into the discharge they inevitably modify 
plasma parameters. In the following we assume that in- 
side the particle cloud the electron temperature remains 
unaffected, while the electron and ion densities are mod- 
ified to keep quasineutrality, Ue + \Q/e\np ~ n^, where 
Q is the particle charge. Furthermore, we assume that 
Tie remains close to the particle- free value no, while n^ 
somewhat increases in response to perturbations from 
the particle component. This approximation is in reason- 
able agreement with numerical simulation results l3, [HI 
regarding plasma parameters inside the particle cloud 
(note that inside the void region in the center of the 
discharge, the plasma density and electron temperature 
can be considerably higher than those in the particle- 
free discharge [l3, [HI ) . Using these assumption we can 
calculate the dependence Q{p^ A) employing the cohision 
enhanced collection (GEO) approximation [19] for the ion 
flux and the OML model for the electron flux to the par- 
ticle. The flux balance condition yields VTi{l + + 
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yi^^ ) = vt, exp(-z), where z = \Q\e/aTe 
is the reduced particle charge, P = (aTe/e^)(np/no) is 
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the scaled ratio of particle-to-plasma densities (the so- 
called Havnes parameter), r = Tg/T^ is the electron- 
to-ion temperature ratio, '^Te(i) = \/T^e{i)l^e{i) is the 
electron (ion) thermal velocity, Aq = Y^Tj/47re^no is the 
unperturbed ion Debye radius, and ii Tn/pciin is the 
ion mean free path {ain ^ 2 x 10~^^ cm^ is the effec- 
tive ion- neutral collisions cross section in argon). The 
interaction potential is assumed to be of Yukawa form, 
the interaction energy between neighboring particles is 
U = (QV^)exp(-A/A), where A = Xo/VTTzP ap- 
proximates the plasma screening length. We neglect long- 
range inverse-power-law corrections to the Yukawa po- 
tential [20I since the ratio = A/ A is not large enough 
in the present experiments (k, 3 and 4 for small and 
big particles, respectively). 

Using the obtained values of A, Q, and A we can es- 
timate variations in U. The results are shown in Fig. 
|4]3. The energy, as a function of p, exhibits maximum 
dit p 11 Pa. The observed trends are in good qualita- 
tive agreement with the results of the structural analysis. 
From the observation that at high pressures the system 
of big particles is near melting (Fig. H^) we can roughly 
estimate the particle temperature, Tp < 5Tn — 0.15 eV, 
by using U/Tp c:^ 7.5 at melting of the Yukawa system 
with = 4: [21]. This estimate is sensitive to the input 
values of Tg and no- To reach equilibrium {Tp Tn) we 
need to decrease Tg or increase no by approximately a 
factor of 7. Unfortunately, Tp is unknown here since the 
particle dynamics is not resolved. Finally, our best esti- 
mate of the relaxation time (defined as a time needed for 
a particle to diffuse over the distance A) yields ~ 1 s (~ 2 
s) for the system of small (big) particles in the fluid state 
near freezing. This is considerably shorter than the char- 
acteristic time associated with external (pressure) per- 
turbations. This implies our complex plasmas quickly 
reach equilibrium static configuration, i.e., they remain 
in quasistatic equilibrium, supporting our interpretation. 

To conclude, let us discuss some specific properties of 
the phase changes observed. To identify the main mech- 
anism responsible for the structural changes let us define 
the relative change of a parameter / between the cases 
of highest (30th scan) and lowest (12th scan) pressures 
as 6f = (/so — /i2)//i2- Normally, one would expect a 
decrease in Q, an increase in no, and therefore a decrease 
in A when p increases. However, the coupling between 
ng, n^, and Up and the corresponding charge reduction 
in dense particle clouds (which is to some extent sim- 
ilar to a reduction of colloidal charge when increasing 
colloidal volume fraction ^]) interferes. As a result, 
in the regime investigated, some increase in A with p 
[6 A 0.11(0.12)] is accompanied by almost constant 
particle charge [Sz 0.01(0.00)] and slight increase in 
A [^A 0.08(0.07)] for the system of smah (big) parti- 
cles, respectively. Thus, an increase in A (i.e. decrease 
in the particle density) is the main factor responsible for 
the melting. On the other hand, reducing pressure com- 



presses the particle system and it freezes. This behavior 
is opposite to the conventional procedure of melting flat 
plasma crystals by reducing the pressure in ground-based 
experiments [3,!23|. The difference is not a consequence 
of the essentially 2D character of crystals investigated on 
Earth, but is rather due to the presence of strong electric 
fields (and, therefore, strong ion flows) required to bal- 
ance the force of gravity. There are effective mechanisms 
of converting the energy associated with ion flows into the 
kinetic energy of the particles. Known scenarios include 
ion-particle two-stream instability , nonreciprocity of 
the interaction due to asymmetric character of the screen- 
ing cloud around the particles ( "plasma wakes" ) 0, [HI 5 
particle charge variations [26|. All these scenarios lead 
to an abrupt increase of the particle kinetic energy at 
pressures below certain threshold value, causing crystal 
melting. The process of melting studied here is quite dif- 
ferent from these plasma-specific mechanisms, but have 
much more in common with those in conventional molec- 
ular and soft matter (e.g. colloidal) systems. 
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